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We investigate the so-called magnification relations of gravitational lensing models. We show 
that multidimensional residue integrals provide a simple explanation for the existence of these 
relations, and an efi'ective method of computation. We illustrate the method with several 
examples, thereby deriving new magnification relations for galaxy lens models and microlensing 
(point mass lensing). 
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I. INTRODUCTION 



T-H ■ Gravitational lensing has proven to be not only important astrophysically, but intriguing niath^matically as well. 

^ \ Mathematical investigations of gravitational lens theory have yielded important results and insights, EraaQEra employing 

^^ 1 techniques and results from such disparate areas as catastrophe theory, differential geometry and Morse theory. 

^^ In this paper, we illustrate how another seemingly unrelated subject, multidimensional residue calculus, applies to 

^j. , gravitational lensing, andpScacifically we explain the origin of certain "magnification relations" that have been discussed 

(•~^ ' in the lensing literature.ElErullj We additionally demonstrate that calculations of these magnification relations are 

f^ , enormously simplified using residue techniques, and illustrate the method by deriving several new results. 
^^ ' This paper is organized as follows. In the remainder of this section, we introduce the relevant terminology of 

gravitational lensing, and describe the magnification relations. In Section O, we express the problem in terms of 

rS ,■ residue calculations at the image positions, and thereby relate it to a residue at infinity. Using one-dimensional 



' ' residue calculus, we derive trace formulas for the magnification relations for a subset of lens models. In Section III, 
^ we consider the general class of lens models, and describe how to perform the necessary multidimensional residue 
integrals. In Section U^, we apply this formalism and derive previously known results, as well as new results. In 
j^ , Section M, we summarize our results and discuss implications. The material discussed in sections El- IV may be 



unfamiliar to astronomers, and so we provide a simple procedure which may be applied to models to obtain their 
magnification relations, without requiring a detailed understanding of the underlying mathematics. 

A Gravitational lensing terminology 

Numerous excellent introductions to gravitational lensing have been written, e.g. Refs.||l|,p|; here we briefly summa- 
rize some of the results and terminology relevant to our discussion. The effects of gravitational lensing can peihaps 
best be understood by considering the time delay of trajectories connecting the lensed source to the observer .Id The 
time delay is a simple sum of two terms, a "geometrical" piece, and a "gravitational" piece. Let 0n be the 3-D 
Newtonian potential, and f/' — {"^1 Dc^) J (/)n dl the projected 2-D potential, where _D is a function of cosmology and 
the source and lens redshifts. Just as V^^n = AnGp, similarly V^ip ex AttGT,, where J^ = J pdl is the surface density. 
In the thin-screen, small deflection limit, the time delay can be written as 

r = roQ|e-/3p-V'), (1) 
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where the normahzation tq depends upon cosmology and the lens and source redshifts, Q is the image location on the 
sky, and f3 is the location the source would havaJiad, had not the lens intervened. Light rays follow null geodesies, 
which can be shown to obey Fermat's principlc.BI Plotting time delay as a function of the two angular coordinates 
on the sky, Fermat's principle demands that images arise at the stationary points of this time delay surface. Setting 
the gradient of the time delay to zero gives the so-called "lens equations" , the (real) solutions of which are the image 
positions. Since the time delay surface can have multiple stationary points, multiple images of a single source can 
arise, and the image multiplicity can depend on the source position relative to the lens. The curves on the source 
plane separating regions of different image multiplicity are called caustics. As a source crosses over a caustic, its image 
multiplicity changes by two, as a pair of images either merge together and annihilate, or are created and move apart. 
The magnification fi, which relates differential area elements of the unlensed source to area elements of the lensed 
images, is simply the inverse of the Jacobian J of the mapping from image coordinates Q to source coordinates f3. 
Since the orientation of an image can be inverted relative to the unlensed source, the magnification can have either 
sign. Much of the Icnsing literature adopts the convention of positive magnifications (i.e. defining /i = 1/|J|); in this 
paper we always take the magnification to be signed, i.e. ^ — 1/ J. 

B Lens models 

Two astrophysically important types of lenses are compact objects (like stars or MACHOs) and galaxies. The 
former class are effectively point masses, and so have lensing potentials proportional to the Green's function for the 
2-D Laplacian, i.e. tp = mlogr where m is proportional to the mass of the point lens. Only the weak-field regime 
is observationally relevant, so the potentials linearly superpose for multiple point masses (as long as they are not 
appreciably separated along the line of sight). Galaxies have more complicated, extended mass distributions with 
correspondingly, complicated lens potentials. In principle, one may decompose the potential into eigenfunctions of the 
2-D LaplacianU 

t/j =2_^{0'mn cos m6 + bmn sin ■m6)r'^. (2) 

This has the advantage that each term in the expansion relates to a corresponding multipolc in the expansion of 
the surface density. In practice, it is necessary to truncate the series due to the limited observational constraints. 
Since galaxies are believed to have roughly "isothermal" p ^ r~^ profiles, such truncated series generally pjnsist of 
variations on the singular isothermal sphere (SIS) potential ip — br. Two examples considered by KochanekEiJ are ths. 
SIS -I- elliptical potential -0 = br{l + 7 cos 26') and the SIS -I- external shear tp — br + ^r^ cos 20. Another variationllZl 

is the singular isothermal ellipse (SIE) ip = bR — b^/siP+]p/q^ with axis ratio q. Other, more elaborate and more 
physically justified models have been employed in the lensing literature; here we will focus on simple models such as 
the above to avoid obfuscating the general method with heavy algebra. In table | we list the models considered in 
this paper. 

Given a model for the lens, the lensed images of a given source may be found by solving the lens equations as 
described above. For simple potentials, the solutions are often analytic. For example, consider the the SIS + external 
shear model. The stationarity equations become 

dr 

— =r- scos{0-es)-b-'jrcos20 = O, (3) 

Idr 

-T^=ssm{e-9,)+^rsm2e = 0, (4) 

where (r, 9) is the image position in polar coordinates, and (s, 9g) is the source position. From this, it is easy to show 
that the quantity u = e*^ satisfies the fourth degree polynomial equation 

^bu + {'-fsv -{ — )u — {sv -\ )u + jb — 0, (5) 

V V 

with V = e*^=. As a quartic equation in u, this can be solved analytically, and from u the image coordinates {r,9) 
follow simply. For many models, it is possible to eliminate all but one variable in a similar fashion and thereby obtain 
analytic solutions, but for more realistic models, the lens equations must be solved numerically. 



C Magnification relations 

As mentioned above, there are in general multiple lensed images for a given source. These images lie at the discrete 
positions satisfying the lens equations for a given source position, and the images have different magnifications and 
orientations. If we restrict ourselves to purely real solutions of the lens equations (i.e. the physically observable 
images) then the number of images changes when the source crosses over a caustic. The number of solutions, of 
course, does not change, but instead a pair of complex solutions-become real (for a source crossing into a caustic) or a 
real pair become complex (for a source crossing out of a caustic). H For certain lens models, such as the simple potentials 
described above, there may exist certain parameter ranges such that all solutions to the lens equations are real. In 
such cases, it has been shown that, there exist interesting and surprising relations among the image positions and 
magnifications. First, Witt & MaoQ considered lensing by a binary microlensing system, involving two point masses, 
and derived the following result : when the image multiplicity is maximized, the sum of the signed magnifications 
of all the images is always 1. That is, X^i M* = 1; independent of quantities such as the lens masses, separation, or 
source position (as long as the source is inside a caustic). This is quite an astonishing result, as the individual image 
magnifications can vary wildly as the source position changes, and even diverge as the source crosses a caustic. Witt 
& Mao derived this result by using resultants to obtain a monic palynomial equation satisfied by ^, and noted that 
the sum of the roots is given by the subleading coefficient. RhieH, using similar reasoning, showed that a similar 
result is true for an arbitrary number of point mass lenses. Dalalu extended this work to galaxy potentials like those 
described above. Again, the magnification relations were obtained by using elimination theory (e.g. Grobner bases) 
to obtain monic polynomial equations in /i; the "total magnification" was then given by the subleading coefficient. 
We summarize these results in-table 1. 

Subsequently, Witt & MaoE3 showed that for a particular class of power-law models, there exist additional mag- 
nification relations, involving not only /i but the image positions as well. They derived this result by separating the 
coupled lens equations into disjoint x and y equations, and then relating /i to the coefficients. Just as previous work 
had derived expressions for J2if^i^ Witt & Mao found expressions for J2if^i^i ^^^^ J2if^iyi- They called these the 
"fc**^ moments", and we adopt their terminology here. For example, they found that the first moment X^iMi^i = 2/3 
for the SIE model. 

The general pattern seen in previous work is that both the total magnification and the higher moments can 
be expressed in terms of the model parameters, and that progressively higher moments have progressively more 
complicated forms involving more of the parameters. The expressions' independence of certain parameters suggests 
some sort of invariant, but clearly not a topological invariant since certain models seemed not to obey any magnification 
relations whatsoever. Indeed, the origin of the magnification relations and their absence in certain models has been 
a mystery. In this paper, we provide an explanation of these relations, and additionally find a method easier than 
elimination theory to derive them. 

D Residue integrals 

As noted above, the lens equations have multiple discrete solutions, both real and complex. While only the real 
solutions have physical meaning, it is instructive to consider the complex solutions as well. In this paper, we will 
henceforth treat the image coordinates Q as complex variables. We are interested in the sum over these discrete 
points of various quantities, such as the signed magnification, or magnification times position, etc. From complex 
analysis, we know that one may relate a sum over discrete points to an integral over a contour encircling those points, 
by choosing an integrand which has poles at those points. For lensed images, which are stationary points of the time 
delay, there is an obvious class of integrands, namely rational functions of the form 

f{x,y) = r. ^ ■ (6) 

OxTOyT 

There are complications, which we discuss later, due to the fact that the integrals here are multidimensional; however, 
the analogy to the one dimensional case should be clear. We need only find the appropriate function g such that / will 
have a residue equaling the quantity we wish to sum over the images, and choose a contour large enough to enclose 
all the images. Now, converting a discrete sum to a contour integral wouldn't seem to be much progress, however 
we can use another idea from one dimensional complex analysis. Recall that by inverting coordinates (mapping the 
origin to infinity and vice versa) one can see that the sum of the residues of poles inside the contour is equal to the 
sum of the residues of poles outside the contour, but with opposite sign. In our case, we are summing over all the 
finite solutions, so the only pole outside the contour is at infinity. This is the essence of the method described in this 
paper: we relate the sum over the images to the behavior of the time delay at infinity, and we simply evaluate the 



residue at the point (s) at infinity. The vahdity of the resulting magnification relations does not depend on the image 
coordinates being real, although of course their physical applicability does. 

In the following section, we further utilize ordinary complex analysis of one variable and derive trace formulas for 



those models for which it is easy to eliminate all but one variable. In Section III, we consider the general case, and 
describe the basics of multivariable residue calculus. In Section |^ we illustrate with specific examples. 

II. TRACE METHODS 

In this paper, we shall focus upon models with polynomial lens equations. Algebraic functions, such as nth roots, 
can be accommodated by introducing an additional variable for each algebraic function along with the polynomial 
equation it satisfies. For example, an equation containing ^Jx^ + 1 is handled by introducing z satisfying z^ — x^ + \. 
The method of this section is useful when all but one variable can be conveniently eliminated from the lens equations, 
e.g. Eq. (pf). Of course, this is always possible in principle, but it may require computer implementation of Grobner 
basis algorithms in practice. Thus, we assume that the x-coordinates of the images are the roots x^ of a polynomial 
equation of degree n, 

n 

f{x)=Y,a,x'^(). (7) 

i=0 

We assume that the signed magnification ^i{x) of an image at x is given by a rational function, 

,-'-^y (8) 

q{x) 

where the denominator has no common roots with f{x). This will necessarily be the case if the coordinates are 
chosen generically, since there will then be at most one image at a given x-coordinate, whose magnification must be 
a single-valued algebraic function of x. We wish to calculate the total magnification 

M = ^/.(x,)=Tr^, (9) 

i 

where the "trace" notation will be explained below. Generically, the roots Xi of f{x) will be distinct; since M is 
determined by continuity when some roots coincide we will always consider the generic case. 

Let A — C[x\ be the ring of polynomials in x with complex coefficients. We call two polynomials g and h equivalent, 
writing g ^ ft,, if they differ by a polynomial multiple of J{x). This sorts the polynomials into equivalence classes, 
and we denote the class containing g by [g] and the set of all equivalence classes by Af . This is of use for our problem 
because all polynomials in a given class take the same values at the Xi and therefore have the same trace. Addition 
and multiplication of classes are well-defined by [g] + [h] — [g + ft], [g][h] — [gh], and Af is itself a ring. 

Now we observe that each class [g] € Af contains a representative which has degree (at most) n — 1. Indeed, the 
relation / ~ implies 



CljiX 



Y,a,x\ (10) 



i=0 



and this can be used to eliminate all terms of degree n or greater from a polynomial g. The resulting polynomial is 
unique, being determined by its values at the n roots Xi. It follows that Af is in fact a vector space of dimension n 
over C. 

Now consider, instead of polynomials, the set R of rational functions of x which are defined at the zeros of /(x), 
and call two elements equivalent if their difference is f{x) times another element. Then all elements of a class have 
the same trace, and the set Rf of equivalence classes is again a ring. In fact, it is isomorphic to Af. An isomorphism 
is obtained by associating to a class [g] in Af the obvious class [g] in Rf. To see that this mapping is invertible, 
we must find a polynomial representative g of an arbitrary class [ft] in Rf. To do so, simply let g{x) be the unique 
polynomial of degree n — 1 satisfying g{xi) = hixi) for all i. Then g — ft is a rational function vanishing at every x^, 
so when expressed in lowest terms its numerator must be a multiple of f{x). Therefore_a e [ft]. 

At this point we can explain the "trace" terminology for the sum over the roots XiX3 The vector space Af has a 
basis consisting of (the classes of) the n polynomials 5i of degree n— 1 defined by 5i{xj) = Sij. Fix an element g € Af 
and consider the linear operator on Af given by multiplication by g. In the given basis, the matrix of this operator 



is diagonal, with entries g{xi). Hence the trace of this matrix, which is independent of the basis and coincides with 
the trace of the operator, is simply Tv g — '^^g{xi). For example, consider Tr x, which just gives the sum of the 
roots. Choosing the basis {1, x,x'^, . . . , x"^^} for Af, the matrix of the operator of multiplication by x has only one 
nonzero diagonal entry, — a„_i/a„, arising from the relation x ■ x"~^ ^ —{an-i/an)x'^~^ + • • • . This recovers the 
stancLaixLi^sult for the sum of the roots of a polynomial and shows how our method generalizes others based on that 
resultClHEle. 

Returning to the problem of computing the total magnification, we see that M = Ti fi = Ti [fj] can be computed 
using any element in its equivalence class. For example, we can choose the unique polynomial representative of degree 
n — 1. Of course, we do not determine this polynomial from its values nixi), since we do not know the Xi explicitly. 
Instead, we seek a polynomial solution of degree n — 1 to the condition (g) defining [//] , 

[n{x)] = [p{x)]. (11) 

This is solved by using Eq. (|o|) to reduce the degree of each side to n — 1, and then equating coefficients. 

To compute the trace we use a formula due to Euler, which|-we derive by means of the residue theorem for complex 
contour integrals. A purely algebraic proof is not difficultE^, but our derivation shows the relevance of residue 
methods and motivates the multivariable generalization which we describe in the following section. Consider the 
contour integral, 

' ^^'(^)^.d., (12) 



27ri/ fix) 

where the contour is a large circle in the complex plane enclosing all the zeros Xi of f{x). The integrand has a pole 
of residue ii{xi) at Xi, and consequently the integral is Tr /i. However, we can also regard the contour as encircling 
the point at infinity and evaluate the integral in terms of the residue there, introducing if we wish the new variable 
u = 1/x to move the point at infinity to the origin. Furthermore, the integral is unchanged if f'{x)ii is replaced by 
any member of its equivalence class. By choosing the polynomial representative of degree n — 1 we need only evaluate 

J-I^dx^^-^^^^^, 0<k<n-l. (13) 

27ri / f{x) a„ 

This proves Euler's formula, 



coefficient in degree n — 1 of /'(a;)/i , , 

coefficient in degree n of f{x) 

where it is understood that the polynomial representative of degree n — 1 is meant in the numerator. This makes it 
clear that the total magnification is determined by the leading behavior of f{x) and f'{x)ii at infinity. 

As an example, we consider a generalization of the SIS + elliptical potential (n = 1 multipole) with an arbitrary 
harmonic, -0 = 6r + -fbr cos m9 . As an aid to clarity, we shall depart from conventional notation, and instead rewrite 
the lens equations so that the variables are x,y, and parameters are denoted a,b,c,.... In this case we define 
x = e'^^,y = r,a = '-f,b = b,c=s,d = e*^= . The lens equations then take the form, 

fix) = mabx^"^ + cd-^x"^+^ - cAx'^-'^ - mab = 0, (15) 

2{y - 6) - c{xd-^ + dx-^) - ab{x"^ + a;"") = 0, (16) 

with the magnification satisfying 

[m,^ab{x"' + x-"") + c{xd-^ + dx-'^)]fi = 2y. (17) 

Eliminating y results in 

[m^abx^"" + cd-^x"'+^ + cdx"'-^ + m^abyi = (18) 

abx^"" + cd-^x"'+^ + 26x" + cdx""'^ + ab. 

Dividing through by x and replacing x^^ in the resulting equation with a polynomial equivalent via the relation 
x~^f{x) ^ produces 



[2m2a6a;2"-i + {m+ l)cd-^T™ - (m - l)cdx'^-'^]^i = (19) 

2a5a;2™-^ + (1 + 771-1)0^-^0;'" + 2bx"'-^ + (1 - 777-i)cdT"-2, 

where the left side is precisely in the Euler form f'{x)ii. Then we immediately have the total magnification invariant 
as 

,, 2a6 2 

M = Ti- n= = — . (20) 

77700 777 

It is no harder to verify that the total magnification is the same for a potential containing an arbitrary finite sum of 
harmonics X^fcLi b^r cos k0; the highest harmonic determines M. Other models are amenable to this method as well. 
In principle one can compute moments by this method, by replacing /i with x'^/i in the relevant equations, but we 
have preferred the residue methods of the next section for moment computations. 

III. RESIDUE METHODS: THEORY 

Let us review the key steps in our contour integral derivation of Euler 's formula in the previous section. First, we 
expressed the total magnification as a complex contour integral of a function having poles at the image locations. 
Second, we converted this to an integral around the point at infinity. This amounts to viewing the complex plane 
as a subset of the Riemann sphere, or complex projective space CP . Our change of variables u = 1/x connects two 
coordinate charts on CP^ centered at the origin and at infinity. Finally, we evaluated the residue at infinity, making 
it clear that the total magnification depends on the behavior of the integrand at infinity. 

These steps all have multivariablc analogs. In fact, there is a well-developed, if little-known, residue theory for 
meromorphic differential forms in several complex variables. This makes it possible to compute directly the total 
magnification and moments for lens models without reducing to one- variable lens equations. The theory is particularly 
effective in the situation of two variables, and we describe it in this case. References include |15ll6Ulj,M. An 
application to chemical reaction rate equations appears in Ref. p9[ , but we are not aware of other physical applications 
in the literature. 



We consider a meromorphic two-form 



g(x,y)dxdy 



Pi{x,y)P2{x,y)' 

on C^ , which we view as a subset of the compact complex projective space CP . Here Pi , P2 are polynomials having 
finitely many common zeros (the image locations) of multiplicity one (as is generically the case), and g is also a 
polynomial. Such a form can be integrated over a 2-cycle, a compact two-dimensional real submanifold of CP . 
Since to is closed (diu — 0), the integral depends only on the homology class of the cycle. For example, in a small 
neighborhood of a common zero of the Pi, we can integrate over the "torus" T : {|Pi| = IP2I = e}, defining the residue 
of uj at this zero: 

Res.= (^)7^.. (22) 

[The standard orientation of T is that specified by the nonvanishing 2-form d(argPi) d(argP2)- That is, T is oriented 
so that dPi dP2 / [IniY P1P2 has positive integral.] We will always denote by J the naive Jacobian of the mapping 
(x,7/) \-> [Pi{x,y),P2{x,y)], namely 



j_ d{PuP2) 
d{x,y) 



d,Pi dyPi 
d^P2 9.P2 



(23) 



This coincides with the physical Jacobian relating corresponding area elements in the source and image planes if 
x,y are rectangular coordinates, and Pi,P2 are the corresponding lens equations, but requires a correction factor 
otherwise. The first key fact we need is that the residue at a nondegenerate zero (one where J does not vanish), 
located say at the origin, is given byta 

Res g^^^^ - A (24) 



which is equal to the magnification of the corresponding image if g is chosen appropriately {g — \ for rectangular 
coordinates). Moments of magnification can be computed by including additional monomial factors in g{x,y). 

Next we need the Global Residue Theoremlla, which states that the sum of all the residues of a meromorphic form, 
such as uj, on any compact manifold, such as CP , vanishes. Note that this sum is over the common zeros of the Pi 
only, so that points where only one polynomial vanishes do not contribute; also the points summed over may depend 
on the choice of the factorization P1P2 of the denominator of uj. The theorem makes it possible to replace the sum 
over the residues at the common zeros in C^ by minus the sum of residues at points at infinity in CP . This is the 
fundamental explanation for the existence of magnification relations in general: compactness relates the sum of finite 
residues to the behavior of the lens equations at infinity, indeed to a finite number of terms in an expansion around 
infinity. It remains to explain how to locate the common zeros at infinity and compute their residues. 

CP^ is conveniently described by homogeneous coordinates [X, F, [/] 7^ [0,0,0], where [XX,\Y,XU] is identified 
with [X, y, U] for all complex A 7^ 0. The points with U ^ can be represented in the form [x,y, 1] and are viewed 
as the subset C^ of finite points. The polynomials Pi{x,y) correspond to homogeneous polynomials Pl^{X,Y,U) = 
U'^'^^^'Pi{X/U,Y/U). Their common zeros at infinity are those having U — 0. These also lie in coordinate charts 
diffeomorphic to C^, described by [l,y,u] or [x, l,u]. In these charts they can be treated just as "finite" zeros are. 
The meromorphic form tu is homogenized so as to have total degree zero, that is, 

^,_ g{X/U,Y/U)diX/U)d{Y/U) 

Pi{X/U,Y/U)P2{X/U,Y/U) ^ °^ 

^ ghudegPi+degP2-dcgg-3^U^XdY - XdUdY - YdXdU) 

ph ph 
^1 ^2 

When n = degg — deg Pi — deg P2 + 3 > 0, the denominator of w'' takes the form V^PiPl^ ■ This creates a subtlety in 
that the above theory applies to a chosen factorization of the denominator into two factors. Thus, we may factor it as 
(C/"P/')(P^) and treat these as the two factors in applying the residue theorem. The common zeros then consist of all 
(finite and infinite) common zeros of the P/', together with any zeros of P2 alone at infinity. The latter might have 
been overlooked in a naive application of the theorem. In the examples we will consider in detail, no such additional 
zeros exist, but we will point out a case where they do. 

Unfortunately, the zeros at infinity are rarely nondegenerate, and their residues cannot be computed using Eq. (p4). 
Instead, they typically lie at singular points of the curves Pi = 0, that is, at least one curve has multiple branches 
meeting at this point.* There is a classical method, dating back to Newton, for finding the branches of an algebraic 
curve P{x,y) = at a singular point, taken to be the origin. The branches are given as Puiseux series, or fractional 
power series, of the form 

00 
j/ = ^a,x"-, (26) 

1=0 

where the exponents ai form an increasing sequence of rational numbers whose denominators eventually stabilize. 
The possibilities for the leading exponent ao are determined by requiring that at least two terms in the polynomial 
P{x,aQX°'°) have the same degree, while the remaining terms have higher degree. Then ao is found by demanding 
the vanishing of the terms of minimal degree. An elegant graphical method for identifying the possible exponents ao 
is provided by the Newton polygon, or diagramOEj. For each monomial x°'y'' appearing in P{x,y), plot the point 
(a, b) in the coordinate plane. Begin at the lowest of the leftmost points (minimize a, then minimize b) and draw a 
polygonal path with vertices at a subset of the points, terminating at the leftmost of the lowest points (minimize 6, 
then minimize a), and choosing each successive segment to have the smallest possible slope (steepest possible negative 
slope). This is the Newton diagram of the polynomial P{x,y). The points lying on any segment of the Newton 
diagram represent terms in P which will have minimal degree if ao is chosen as the negative reciprocal of the slope of 
that segment. The next term in the Puiseux series can be found by applying the same procedure to P{x, a^x""^ + y), 
and so on. One term is often enough for computing the total magnification; higher moments require more terms 
in general. The denominators of the exponents stabilize at the stage where the Newton diagram has only a single 
segment. The curve P{x, y) — has at least as many branches at the origin as there are segments in the Newton 
diagram, and has more if there are multiple solutions to the equations for the coefhcients ai . 

Consider one particular branch X of the curve Pi{x,y) = at a singular point, given by a Puiseux series y = 
aoa;™'" + • • • , with m, n relatively prime. On this branch draw a small circle C around the origin; its projection on 



*Notc that these curves have one complex, but two real dimensions, and so are real surfaces. 
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the x-plane must wind n times around the origin. Now construct a 2-torus 6C by "thickening" C: at any point p of 
C take a plane transverse to X and a small circle in this plane with center p. As p moves around C this circle sweeps 
out the torus SC. We may construct such a torus for each branch of the curve. These are called Leray tori, and the 
"thickening" operator S is the Leray coboundary. 

Our objective is to compute the residue integral L, uj at the origin. We can work entirely within a small ball B 
around the origin. The residue depends on the homology class of T in H2{B — {P1P2 =0}). As we explain in the 
Appendix, this class is, up to sign, the sum of the classes of the Leray tori constructed on the branches of either of 
the curves Pi = or P2 = 0. Therefore the residue is the sum of the integrals of ui over either set of Leray tori, with 
appropriate orientation. 

The integral over a Leray torus SC lying on a branch of Pi (a;, y) — given by a Puiseux series y — p{x) is computed 
using the Leray residue formula. We give this formula under the assumption that dyPi does not vanish on the given 
branch, which amounts to assuming that y ~ p{x) is a branch of an irreducible factor of Pi which appears to the first 
power only (the analog of a simple rather than multiple pole). The general formula can be found in Ref. pq |. Our 
case reads 



1 f g dx dy 

27ri 



SC 



P1P2 



g dx 



P2 dyPl 



y=p{x) 



(27) 



which is proved as follows. To integrate over SC, we can integrate first over the circles in the planes transverse to C, 
then over C itself. Since these circles are centered at Pi = 0, we can change variables from x, y to x, Pi by means of 
dPi dx = —dx dy dyPi and then integrate over Pi by means of the one-variable residue theorem. The evaluation of 
the residue at Pi = by substituting y — p{x) leaves another one- variable residue integral to be performed. For this 
one must keep in mind that the cycle C may wind around the origin several times in the x-plane. 

The Leray residue formula holds with obvious notational and possible sign changes if the Leray torus lies on a 
branch of P2{x, y) — 0, or if we choose to eliminate x rather than y in favor of a P^. 



IV. RESIDUE METHODS: EXAMPLES 



In this section we apply the residue methods just developed to various lens models. The first is a generalized n = 2 
multipole model with potential ip — br + {^/2)r'^ cosm0. We again redefine notation to clarify the computation, 
writing the variables as x, y and the parameters as a,b,c, . . . . Let x = e^^,y = r,a = j,b = b,c = s,d = e*^= . The 
lens equations then take the form. 

Pi == cd-ix"+i - cdx™-i + y ay(a;2™ - 1) = 0, (28) 

P2 = 2(y - 6)x" - cd-^x^+i - cda;"-i - ay{x^"' + 1) == 0. (29) 

Because x, y are not rectangular coordinates, there is an extra Jacobian factor in the magnification, and we find 

„2m-l 



/" = 



4yx 



J 



Consequently the total magnification is given by the residue sum, 

4yx^"^~'^dxdy 

images 



M= ^ Res 



P1P2 



(30) 



(31) 



The residue theorem relates this to the sum of the residues at points at infinity in CP , which are found from the 
homogeneous forms of the Pi : 



cd-^X"'+^U"' - cdX"'-^U"'+^ 



-aYiX^' 



C/2™), 



P^ = 2A"(rt/'" - 6[/™+i) - cd-^X"'+^U"' - cdX"'-^U"'+'^ - aY{X^"' + f/^'"). 



(32) 
(33) 



The common roots at infinity are those with [/ = 0, and there are two: [X, Y, U] = [1, 0, 0], and [0, 1, 0]. The total 
magnification is minus the sum of the residues at these points of 



4y^2m-l[/2-m+2 rf(X/C/) d{Y/U) 
^1 ^2 



(34) 



• 1 >- u 

m m + 2 2m 

FIG. 1: The Newton diagram for P^^^ in the n = 2 multipole model The degree of a monomial in y (u) is plotted 
vertically (horizontally). The case ttt, = 3 is shown. 

Consider first the point [1, 0, 0], which we examine in the affine chart [X, Y, U] — [1, y, u], where 



P^=' = crf-iy™ - cdw™+2 + -ay(l - u^m)^ 

P2^=i = 2M™(y - bu) - cd-^u"" - cdu"'+'^ - ay{l + u^™) 



(35) 
(36) 



and we need the residue at the origin of 4j/u^™~-'^ du dy / P^^^^ P^^^ . 

The Newton diagram for Pi is shown in Figure 1; there is a single branch on which to leading order y '^ w™, as is 
easily verified by solving P^^^ — for y. The Leray formula evaluates the residue as the one-variable residue of 



4:yu'^"'-'^du 

pX=la pX=l ' 
^2 ^y^l 



(37) 



where y ~ u™. But it is easily seen that the leading behavior of this 1-form near m = is u^™~^ du, so that there is 
no pole and no residue for tti > 0. 

We examine the remaining root [0, 1, 0] in the chart [X, Y, U] = [x, 1, u] where 



P^=^ = cd-^x^+iu" - cdx™~iu™+2 + -a(x2™ - u^™)^ 



jY=l 



2x'"u™(l - bu) - cd-'x'^+'u"' - cdx"'-'u"'+' - aix^"" + u^™) 



and we compute 



Res 



4a;2™-iu2m-i 



dx du 



pY=l pY 
^1 ^2 



(38) 
(39) 

(40) 



The Newton diagram for P^ is shown in Figure 2, and gives the leading behavior u ~ /x+- • • . With this behavior, 
the lowest-order terms of P^^^ vanish iff P™ = 1, so that there are 2m branches with Ip = exp{ip7T/m), p — 1, . . . , 2m. 
With the Leray formula, the total magnification becomes 

4^2m-ly2m-l ^^ 



Y^ -Res 

U — lpX-\ 



Pf=i9„Pf=i 



(41) 



The leading terms are readily identified, and indeed the behavior is as dx/x, with residue 



2m 



2 ^ 1_ 



^-^ m^ail™ — a) m?a ^-^ (—1)^ — a m{\ — a?) 
p— 1 ^ ^ ' p—i ^ ^ ^ ^ 

In terms of the original parameters of the model, 

4 



M 



to(1 — 7^) ' 
an attractive generalization of the known resulto for m = 2. 



(42) 



(43) 
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FIG. 2: The Newton diagram for P-[ in the n = 2 muhipolc model. 



Once the branches of the Pi have been identified, it is easy to modify the calculation to compute moments rather 
than total magnification. For example, let us compute the first x-moment, Tr x/x. (We continue to use the Tr 
notation for the sum over the images.) Since x = e'^ in terms of the physical variables of this model, from the real 
and imaginary parts of the result we can obtain the moments weighted by cos 9 and sin0, in the case that all images 
are real. In homogeneous coordinates this simply gives an additional factor X/U in the residue in Eq. (pi]). At the 
point [1, 0, 0] this adds a factor 1/u to the one-variable residue (37), changing the behavior to u^™"^ du. There is still 
no contribution for ra>\. At the point [0, 1,0], there is an extra factor x/u — l^'^ + ■ ■ ■ in the residue (^l|), changing 
the contribution to 



E' 



I (JjiipXi jj 



2m 

■E 



-ipTT /rn 



^ (-i)p 



(44) 



which can be evaluated in closed form if desired. 

From the real part of Tr yx[i we can obtain the moment of r cos d, the physical x-coordinate of the image. This 
leads to a factor YXjU'^, which worsens the singularity at [0, 1,0] to a double pole, requiring an additional term in 
the Puiseux expansion to obtain the residue. The result, for the true external shear model ttt. = 2, is 



Tr [IT cos d — 



ZXt; 



(l-7)(l-72)- 



(45) 



This example illustrates the general situation. Higher moments produce extra monomial factors in the residue 
expression. In general this will worsen the singular behavior at points at infinity, although this may not occur for 
certain branches, such as those at [0, 1,0] in the example of Tr xfx. This will have two effects: points which do not 
contribute to the total magnification generally will contribute to higher moments, and more terms in the Puiseux 
expansions will be required for higher moments. Both effects will result in higher moments being given by more 
complex expressions, with more model parameters contributing. nn 

We turn to a second example, the Singular Isothermal Ellipse (SIE) potentiallj'EBl. For this model ip = bR = 
h\/ x^ + y^q'^ ^ where x, y are rectangular coordinates and 6, q are parameters. The lens equations are 



bx 

'^x ^ X Xs ~^ — U, 
rC 



y-Vs 



by 

q^R 



(46) 
(47) 



and the magnification is given by 



M-^ = 



Txy Tyy 



(48) 



11 



Algebraic manipulation leads to the polynomial equations, 

pi == q'^XTy - yr^ = q'^x{y - y^) - y{x - Xs) = 0, (49) 

P2 = R^tI - 2bRxT^ = [x - Xsf{x^ + ?/g"^) - b^x^ = 0. (50) 

However, these equations have the extraneous solution (x, y) — (0, 0), which does not satisfy the original lens equations. 
This can be eliminated by substituting y = wx^ and adopting the modified equations 

Pi = pi/q X ^ wx — ys — rjwx + rjwXs, (51) 

P2=P2/x^=b^-{x-Xsfil + rjw^), (52) 

where we have set t] = q^^ . Relating the naive Jacobian J of the Pi with respect to a;, w to the Hessian of r gives the 
magnification in the new variables, 

2x{x - Xs){l + W) , . 

fj,= J . (53) 



The moments of the magnification with respect to x are given by 



Tr a; ^ = 2^ Res — — . (54) 



In homogeneous coordinates we have 

Pi = (l - 'n)XW - ysU"^ + riXsWU, (55) 

P^ = b^U'^-{X-XsUf{U^+'qW^), (56) 

and we need to compute 

Tr ,x^ . ^ Res 2X^^^i7^-^(X - xM)iU^ + ,W^) djX/U) djW/U) ^^^^ 

There are two common zeros of the P/* at infinity, namely [X, W,U] — [1,0,0], [0,1,0], and in the corresponding affine 
charts we have 

Pf^^ = (1 - '>l)w - ysu"^ + Tjxswu, (58) 

Pf =1 = (1 - r])x - ysu^ + TlXsU, (60) 

P^-^ ^b'^u'^^{x-xsuf{u^ + ri). (61) 

The residues at these zeros may be computed via the Leray formula applied to the branches of either Pi or P2 ; we 
have done both computations and the latter seems slightly simpler. In each case the branches can be determined 
directly without appealing to the Newton diagrams. At [0, 1,0], the equation P]^""^ = is solved by 

V? 

x:^XsU±bri-^/'^u^{l + —)-^/'^, (62) 

giving one branch for each choice of sign. The contribution to Tr jix*' from one branch is 

2x^^^{x — Xsu){v?' + rj) dxdu 

^^^ yk+lpW=lpW=l • ("3) 



Applying the Leray formula gives 



which simplifies to 



2x''+Hx-Xsu){u^ + iq)du 

^k+lpW=lQ^pW=l ' ^o^i 



u'^^^ du 
^[x-ysU^ -ri{x- Xgu)] 



T^ X au u ,„_. 

R-es .+ir ....2 „f^ ^ .,M = ^s , (65) 
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where only the first term in the series expansion of x was required. The two branches at this point thus contribute 
2x^ to the moment of x^ . 

At the remaining point [1, 0, 0], P^^^ = is solved by 




± — (1 --&'«' + •••), (66) 



{l-XsuY ^' 2 



and we need 



2{1 - Xsu){u'^ + r]w'^)dudw 
I^es ^k+ipx=ipx=i ' (67) 

summed over the two branches. Integrating over w using the Leray residue formula and simplifying yields 

{u^ + rjw'^) du 

-Res , —r- 7 —, (68) 

rywu'^+^^l — XsUj[[l — rjjw + rjXsWU — ys" J 

where u^ + rjw^ — iP'u^ + • • • . Since w ^ u, there is no pole for k < 2: the total magnification and first moment are 
simply given by the contributions from [0, 1,0]. There is an additional contribution to the second moment given by 
-2b'^/{l-r]^) =2b'^q^/{l~ q'^), which agrees mth the result of Ref. ||l^. 

We have verified the results of Witt & MaoE£l through the third moment. Those authors noted that the moments 
of X were independent of j/s to this order, and that "this seems remarkable" in view of the dependence of /i on this 
parameter. The explanation is that the term ysU^ in P^^^ does not contribute to the residue for low moments; indeed 
it does contribute to the residues for the third moment but its contribution cancels between the two branches. The 
fourth and higher moments do depend on t/s. 

Lastly, we note that it is entirely straightforward, and no more work, to generalize the calculation to include an 
arbitrarily oriented external shear term. The first moment, for example, takes the form 

2 2 

Tr ^x= {xs + -liXs + 722/s), Tr ^ly = {y^ + 722;^ - 71^^) (69) 

where 71 = 7 cos 26*.^, 72 = 7 sin 20^, and d^ is the orientation angle of the shear (see Table ||). 

As our final example we consider microlensing due to a collection of N coplanar point masses. We adopt Witt's 
complex notation, writing z = x + iy for the position of an image and w for the position of the source. The lenses 
have masses rrii and positions Zi. The lens equations are 

z-u;-f:^^=0, (70) 

and its complex conjugate; when we complexify the coordinates x, y of an image, z and z become independent variables 
and the conjugate equation becomes an independent condition as well. The observable (real) images are those for 
which z is the conjugate of z. Clearing denominators, we set 

Pi - (z - w) Yliz -z,)-Y, m, [](z - z,), (71) 

i i j^i 

P2 = (Z - W) 1[{Z -Z,)-J2 ^^ n(^ - ^^■)- (72) 



For the magnification we find 



lJi = J-^\{{z-z,){z-z,), (73) 



where J = d{Pi, P2)/d{z, z). For the /cth moment of magnification, wc must compute 



fe_-D„„ z'^Yliiz- Zi){z- Zi)dzdz 



Tr ^lz'' = Res — ^^^ — . (74) 

P1P2 
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N >' 



N + 1 <> 




FIG. 3: The Newton diagrams for P^ (left) and P2 (right) for N point masses. The case iV = 4 is shown. 



The homogeneous polynomials are 

Pi''- = (Z - wU) YliZ - z,U) - J2 ^^U^ X{{Z - z,U), 



j^i 



P^ = {Z~ wU) Y[iZ - z,U) - Y, m^U^ Y[iZ - ZjU). 

i i j^i 



(75) 
(76) 



Setting f7 = 0, we obtain Pf = ZZ , Pj = ZZ^ so there are two common zeros at infinity, [Z,Z,U] = 
[1,0,0], [0, 1,0]. However, we can say more: each of these zeros has multiplicity N. Since the homogeneous polyno- 
mials each havepdegree A^ + 1, the number of finite common zeros will be (A^ + l)'^ — 2N ~ N^ + 1, in agreement with 
previous resultsB. 

Dehomogenizing the polynomials at these points, we find 



Pi = (z — Wu) 11(1^ UZi) — U 2_, "T-J I I (1 ^ "^Zj)^ 



3^1 



P2 = (1 — wu) \\{z UZi) — U 2_^^^i I I (•^ ~ ^-^j)] 



j¥« 



Pi = (1 — wu) I I (^ — UZi) ~ U 2_^^^i I I (•^ ~ "^Zj), 



]^l 



(77) 
(78) 
(79) 
(80) 



P2 = (z — Wm) 11(1 — UZi) — M >^77T,j I I (1 — UZj). 

i i j^Li 

Because of the complex conjugation symmetry of these expressions, it suffices to examine the branches of, say, P^ 
deduce the others. The Newton diagrams of these are shown in Figure 3. In each case, the Puiseux series are ordinary 
power series. For Pi there is a single branch u = z/w+- ■ ■ , while for P2 there are N branches u = z/ Zi—{mi/ zf)z'^+- ■ ■ . 
It is now straightforward to compute the residues 



z=i 



to 



Res 



Zku2-k jji^z _ ZiU){Z - z^U) d{Z/U) d{Z/U) 



r>h TDh 
^1 ^2 



(81) 



and determine the moments. The first two are Af = Tr /i = 1, and Tr z/i — 'J^+X^i -'"-. ■ The single branch contributes 
w^ to the fcth moment, while the contribution of the N branches becomes progressively more complicated. 

We have also considered a generalization of the model to include external shear. This brings the lens equations to 
the formO 



z — w — 7Z - 



N 

E 



= 



Z- Zi 



(82) 
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Model xf) Y.i^J■i Y.jfJ'iZi 

point masses J^^ m; log |9 - 6; | 1 Zs + J2i IJ^ 

point masses + shear tppm + -^r^ cos 26 1/(1 — 7^) n- '^y^ (-^° + 7^^ 
SIE bR 2 ^ 2zs 



^3 



SIS + elliptical hr + 'ybr cos 20 1 Zs + 272^ — „„.| ^ 

SIS + shear 6r+^r2cos26» 2/(1 -7^) ^_1^(^, +7^) 

SIE + shear bfl + f r^ cos2(6) - 6*^) 2/(1-7^) ^^ii^j (z« + 76^'^''^,) 

TABLE I: Model lens potentials and results. The t^^ lens has mass mi and position Qi. For the galaxy potentials, 
with only one lens, we choose coordinates centered on the lens and oriented along the ellipticity or shear axes. Here, 
the image position Q = {xiV) — {r cos 9, r sin 9) and source position /3 = {xs^Vs) — {scos9s,ssu\9s)- We use the 
complex notationci z — x + iy, and similarly for Zs,zi. Also, R = \J x'^ + ip j (p' is an elliptical radial coordinate, 7 is 
the strength of the shear, and q is the axis ratio. 

and its complex conjugate. Defining P^ as before, their zeros at infinity are [1, 0, 0] and [7, 1, 0] for P/', and [0, 1, 0] and 
[1,7,0] for P2 ■ For 7 ^ 0, 1 there are no common zeros at infinity. Factoring the denominator of w'' as {U"Pi){P2), 
the residue theorem gives the sum over the finite common roots as minus the sum of residues at the infinite roots of 
P^. We find, for example, the total magnification 1/(1 — 7^), and the first moment 



w + ^w 



Tr^M^T^ -^- (83) 



Note that the first moment's dependence upon the lens' positions has vanished. 

V. DISCUSSION 

We have introduced in this paper a new framework for analyzing gravitational lens models. The use of residue 
integrals makes clear the origin of the magnification relations, and facilitates their computation for a wide class 
of model potentials. We have also applied this method to a series of models, confirmed and extended previous 
results, and provided new magnification relations for several models. Although multidimensional residue calculus 
may be unfamiliar to readers from astronomy, one may follow a simple procedure to perform the necessary integrals. 
Basically, the procedure is as follows : 

1. From the stationarity equations, construct two polynonfials Pi,P2 that simultaneously vanish at (and only at) 
the image positions. 

2. Define the "Jacobian" Jp = dct[d{Pi,P2)/d{x,y)], and define g{x,y) — Jp{x,y)fi{x,y) where n is the magnifi- 
cation. 

3. Change to homogeneous coordinates: {x,y) -^ {X,Y,U) with x = X/U,y = Y/U and homogenize the polyno- 
mials by multiplying each by the factor f/'^'^s ^. Also, multiply g by U'^°s ^i+deg P2^ 

4. From the homogenized g, Pi, P2, construct the 2-form {g / PiP2)d{^)d{jj-) , and factor the denominator into two 
groups; usually, the grouping Pi,P2 suffices. (Henceforth Pi,P2 shall refer to the two groups, not the original 
polynomials.) If the denominator contains explicit factors of U then redefine one of the polynomials, say P2, to 
contain these additional factors. 

5. For U — 0, find the points {X,Y) where Pi,P2 simultaneously vanish; these are the roots at infinity. 

6. Pick one of the polynomials - say Pi - and determine the behavior of Pi = in the vicinity of each common 
root. First define coordinates for the neighborhood of the root. For example, if the root is Y = 0,U — 0, then 
a good choice would be {X,Y,U) — {l,y,u). As discussed, there are in general multiple branches of Pi =0 
meeting at the root at infinity, each parametrized by a Puiseux series y = ^^ atu"* . The ai's can be determined 
from the Newton diagram, and by substituting in the specified a^'s one may solve for the coefficients a^. Usually, 
only the first one or two terms in the series are necessary. 
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7. Construct the quantity 



g{u, aiu"^ + a2U° 



^{u, aiM"i + ...)P2(u,aiu"i + •••) 



and pick out the term ^ u ^ . The coefficient of u ^ is the contribution for this branch; summing over all the 
branches gives the residue for each root. 

8. Repeat this procedure for all the roots at infinity, and sum their residues. Negating this quantity gives the sum 
of the residues at finite poles (the images). 



The above is of course just a rough outline; section III describes the method in full detail. Following this procedure, 
the results listed in Table | can be reproduced with minimal effort. As mentioned earlier, the Newton diagram method 
of determining the branches near each pole is effective only if the lens equations can be brought into polynomial 
form, such as for these simple potentials. For more complicated models, introduction of additional variables may 
be necessary to handle fractional powers. Our method may not be applicable to models involving transcendental 
functioris, in particular functions with essential singularities at infinity. 

Dalala and Witt & Mac£3 have considered the applicability of such magnification relations to real gravitational 
lenses. For galaxies, DalalS has shown that these relations can be an aid in fitting models to lensed objects, or can be 
used to rule out models a priori. Witt & MaoE3 have shown, however, that reliance upon simple galaxy models can 
be misleading, when applied to realistic galaxy potentials. This limits the applicability of magnification relations to 
making statements about models, as opposed to statements about the lenses themselves. For microlensing, however, 
there is no doubt about the accuracy of the point mass approximation, and as such our derived magnification relations 
may be considered exact. Witt & MaoQ have already shown how the total signed magnification ("zeroth moment") 
can be used for binary microlensing to set lower limits on the overall unsigned magnification, useful for example for 
detecting source blending. Although the multiple images of a microlensing event cannot as yet be resolved, precluding 
the present-day experimental verification of our prediction regarding the first moment, we are hopeful that the future 
advent of space-based interferometers will allow our microlensing formulae to be tested observationally. 
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APPENDIX : THE HOMOLOGY CLASS OF THE TORUS T 

Here we explain why the torus T defining the local residue at a singular point, taken to be the origin, is homologous 
to the sum of the Leray tori constructed from the branches of either of the polynomials Pi{x, y), and how to determine 
the correct orientations. We |ase grateful to Eduard Looijenga and Peter Teichner for explaining the topology to us. 
See also section 2.2 of DimcacJ. 

We are working locally in a closed ball B around the origin, and we denote by X the zero locus P1P2 ~ within 
B. The various tori define homology classes in i?2(-B — X). It is known tha|t-|X,is topologically a cone with vertex 
the origin and base X n dB, which is a linked collection of topological circlegl3'El. Denote the several branches of X 
as Xi, and the Leray torus built on a given branch as 6Ci. Let aij be a path running from the origin to dB along Xi, 
and returning to the origin along branch Xj. ApSubset of these paths forms a basis for the relative homology group 
Hi{X,dB). Furthermore, by Alexander dualityllj, an element of H2{B — X) is uniquely determined by its linking 
numbers with these paths (indeed, with those in a basis alone). (The linking number l{c^,c^) of a 1-cycle with a 
2-cycle in a 4-ball is the intersection number of c^ with any 3-manifold having boundary c^.) 

The Leray tori have linking numbers 

l{aij,6Ck) = Sik - Sjk- (Al) 

Indeed, the intersection of aij with a solid torus bounded by 6Ck is the intersection of aij with Ck, which is one point 
if the outward segment of aij lies on branch Xk , and one point (with opposite orientation) if the returning segment 
does. 
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The sum J^p ^^k of Leray tori built on branches of Pi =0 therefore has hnking number with aij equal to +1 if Xi 
is a branch of Pi and Xj is a branch of P2j —1 if vice- versa, and if Xi and Xj are branches of the same polynomial. 
The sum "Y^p SCk has the negatives of these linking numbers and therefore represents the same homology class but 
with opposite orientation. 

It remains to show that T has the linking numbers of J^p ^Ct- A solid torus bounded by T is given by {|Pi| < 
e, IP2I = e}- This meets any branch Xk of Pi = given by a Puiseux series y — p{x) in the locus \P2{x,p{x))\ — e 
on Xk- This set is topologically a circle around the origin, having intersection number +1 with a path radially 
outward from the origin. The solid torus meets no branch of P2 — 0. Therefore, radially outward (inward) paths on 
any branch of Pi contribute +1 (—1) to intersection numbers with this solid torus, while paths on branches of P2 
contribute nothing. This duplicates the linking numbers of ^p SCk- 



The orientation for T used in this argument is indeed the standard one prescribed in section III. Using the Leray 
residue formula to evaluate 

2m) k P1P2 ^ ' 

produces a positive contribution from every branch of Pi . 
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